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, We describe and characterize rigorously the chaotic behavior of the sine- 

Gordon equation. The existence of invariant manifolds and the persistence of 
homoclinic orbits for a perturbed sine-Gordon equation are established. We ap- 
ply a geometric method based on Mel'nikov's analysis to derive conditions for 
the transversal intersection of invariant manifolds of a hyperbolic point of the 
""Q ' perturbed Poincare map. 
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1 Introduction 

Homoclinic orbits have long been identified as a possible source of chaos in nonlinear 
finite-dimensional dynamical systems. For example, when an integrable 2-dimensional 
system, containing a smooth connection between the invariant manifolds of a saddle 
fixed point is perturbed, this connection is generally broken and the invariant manifolds 
intersect transversally. This yields isolated homoclinic orbits of the perturbed system, 
near which one expects to find shift map embeddings and chaos. 

In the theory of finite-dimensional dynamical systems, a standard method to convert 
numerical experiments of chaotic dynamics into a rigourous mathematical description 
is to construct Smale horseshoes in a neighborhood of a homoclinic orbit [25] and then 
use these horseshoes to establish the existence of an invariant set on which the motion is 
topologically equivalent to a Bernoulli shift on a number of symbols. Such constructions 
begin with the existence, persistence and breaking of homoclinic orbits, which in turn 
are obtained through Mel'nikov analysis. 

However, in the case of infinite dimensional dynamical systems the situation is more 
complicated [24]. As is well-known, realistic physical systems are usually modelled by 
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nonlinear PDEs, whose chaotic behaviour is generally very difficult to study analytically. 
In the last 10 years, D.W. McLaughlin and collaborators have studied the periodic spec- 
tral transform of integrable soliton systems and constructed representations of global 
objects such as homoclinic orbits and whiskered tori, extending the Mel'nikov analysis to 
the NLS equation [15], [17]- [19], as well as a system of coupled NLS [22], [23]. Recently, 
G. Haller has been developing an alternative prespective, under certain conditions, he 
proved the existence of manifolds of multi-pulse Silnikov type orbits homoclinic to the 
critical torus of NLS equation. He has applied to finite dimensional discretizations of 
the perturbed NLS equation [8], [9] and to the PDE case [10], [11]. 

In this article we describe and characterize with mathematical precision the chaotic 
behavior of the perturbed sine-Gordon equation (SG) 

Casel u t t — c 2 u xx — sin-u = e{gi{u) — au t ) (1.1a) 

Casell u tt — c 2 u xx — sin-u = £(<?2(0 — au t ) (1-lb) 

with even periodic spatially boundary conditions, e a small perturbation parameter, 
< £ < 1. The perturbation term g\{u) is a smooth nonlinear function satisfying 
<7i(0) = 0, d u gi(0) = and g 2 is time-periodic with period T. 
The unperturbed sine-Gordon equation 

u tt — c 2 u xx — sin-u = (1.2) 

is an infinite-dimensional Hamiltonian system [24]. We denote by H 1 the Sobolev space 
of functions of x which are L-periodic, even and square integrable with square integrable 
first derivative on [0, L) 

H 1 = | u = (-u, u t ) : u(— x) = u(x) = u(x + L), J \d x u\ 2 dx < oo j = X (1.3) 

The system of sine-Gordon models the phase difference between two superconducting 
layers in a Josephson junction, [21]. Equation (1.2) can also be thought of as a model 
for the continuum limit of a chain of coupled pendula, including damping and driving, 
with u(x,t) the angle from the vertical down position of the pendulum at the point x 
along the chain. 

On the phase space H 1 , the Hamiltonian of the equation (1.2) is given by 
H : H 1 — >R, H(u) = £ {^{u 2 t + c 2 u 2 x ) + (1 + cosw)) dx 

with an infinite number of commuting constants of motion. The common level sets of 
these constants of motion are generically infinite dimensional tori of maximal dimension. 
These constants of the motion have linearly independent gradients, at certain critical 
tori, which become linearly dependent when the critical tori have dimension lower than 
maximal (but otherwise arbitrary). Critical tori can be either stable or unstable. In the 
unstable case, the constant level sets have a saddle structure in a neighborhood of the 
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critical torus. An unstable critical torus has an unstable manifold of infinite dimension. 
If, a phase point on this manifold approaches the critical torus as t — > ±00, such a 
phase point is said to lie on an orbit which is homoclinic to the critical torus. 

Bishop et al. [1] have studied the sine-Gordon attractor numerically for g2(t) = 
Tcosut and the parameter region 0.1 < eT,ea < 1. In that region the so-called 
"breather" solutions (i.e localised in space and periodic in time) interact and quali- 
tatively the dynamics appears to be well described by a four mode truncation [8], [14]. 

Ercolani et al., [2] have explored the existence of homoclinic orbits for the unper- 
turbed system. To this end, they used the fact that system (1.2) admits a Lax pair rep- 
resentation and derived explicit formulas for all homoclinic orbits in terms of Backlund 
transformations. In these formulas, the use of the imaginary double points of the spec- 
trum of Lax operator is very important, since it describes the homoclinic orbits in 
infinite-dimensional space. 

However, it is very important to study the behavior of the trajectories of the per- 
turbed problem (1.1a), (1.1b). Especially, it is interesting to ask: under which pertur- 
bations do the homoclinic orbits persist and how does one establish their persistence? 
In the finite-dimensional system, the method which is most often used to answer such 
questions is the one established by Mel'nikov [20], [25]. 

In this paper, we extend and apply this geometric method for the perturbed system 
(1.1a), (1.1b), to establish the persistence of homoclinic orbits through the transversal 
intersection of invariant manifolds of singular fixed point. In section 2, we describe the 
symplectic structure and briefly outline the theory of the e = (integrable) case. In fact, 
in the integrable setting, the associated spectral problem has been used to classify all 
instabilities and to construct representations of their homoclinic orbits and we derive an 
analytic expression for the gradient of the Floquet discriminant. Section 3 is devoted to 
describe the invariant manifolds for the perturbed system (1.1a), (1.1b). We study the 
spectrum of the linearised system in the neighborhood of the uniform solution u — (0,0) 
lies on the circle S. The linearized operator admits two pair of real eigenvalues Aq , Af 
and infinite number of imaginary eigenvalues A^ = ±iaj, j>2, thus the uniform solution 
(0, 0) of the sine-Gordon equation in the whole space becomes a singular point of type 
saddle-focus. We rewrite our system in a more convenient form and prove an invariant 
manifold theorem for the perturbed sine-Gordon equation (1.1a). Also, we construct the 
corresponding Poincare map for the time-periodic perturbed system (1.1b) and analyze 
the dynamics on the perturbed manifolds. 

We derive, in section 4, a formula for the distance between the codimension two 
center-stable and unstable manifold of the critical torus S of the perturbed system. 
We prove the existence of persistent homoclinic orbit to the saddle-focus point £ = 
+ 0(e) E S for the system (1.1a) and for (1.1b) with g 2 (t) = Tcosut, this problem has 
been studied numerically by Bishop et al. [1]. The derivation of the distance between 
the manifolds is an extension of Mel'nikov methods in infinite dimensional space. We 
prove the convergence of the Mel'nikov integrals in Case I and II. The existence of simple 
zeros of the distance function implies a transversal intersection of those manifolds, as 
long as certain additional non-degeneracy conditions are satisfied. These simple zeros 
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correspond to orbits that are asymptotic to the saddle-focus point £ in backward time 
following the curve which lies in the unstable manifold and in forward time following 
the curves on the center-stable manifold. 



2 Discussion of the Integrable System 



Consider the nonlinear wave equation (1.2) on a segment of length L with even periodic 
boundary conditions in a symplectic space X (1.3) with scalar product ( . , . ), endowed 
by a weakly symplectic form fl given by 



f L 

Q(u, u 1 ) = / (ulu — u t u l ) dx, Vu, u 1 E X 
Jo 



(2.1) 



Let us define a domain V C X and a function / G C l (V) and let V Q / G X be the 
gradient of / with respect to the inner product in X 



1 r L 

Df(a)w := — / V a f(a(x))w(x)dx 
L Jo 



Denoting by 



H(u) = £ {^{u 2 t + c 2 u 2 x ) + (1 + cosw)) dx G C\X) 



the Hamiltonian vector field Vh with Hamiltonian H and the map Vh '■ V 
that Vh = JVif(u), the equation (1.2) takes the form 



X such 



/ 1\ 
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\ c u xx + sin u J 



(2.2) 



Defining the Poisson bracket 



r L /5F _ 5F_ 5G\ 

' jo \ 5u Su t 8u t 5u J 

we note that the evolution of any functional F, under the SG flow, is governed by 

dF 

{F,H} (2.4) 



dt 

Obviously, the Hamiltonian H is conserved by the SG flow. The SG equation admits an 
infinite family of conserved functionals in involution with respect to the Poisson bracket 
(2.3). This fact allows the SG to be solved with the inverse scattering transform [4]. 

Let us denote by D(V H ) = {u e V, Vff(u) G X} the domain of V H in V. We give 
the following definition which is useful for the study of invariant manifolds, contained 
in section 3. 
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Definition 2.1 Let V\ C V be such that for every Uo € T\ there exists a unique solution 
u(t) = $*(u ), t 6l of (2.2) with initial condition u = u(0). The set of the mappings 

&:Pi — >V, u — >$'(u ) tel 

is called the flow of equation (2.2). 

We are interested in studying the dynamics of the perturbed system (1.2) on a 
submanifold, X, of codimension 2, in the infinite-dimensional phase space X, given by 
two functionally independent equations 

X = {u E X : u x (x,t) = 0} (2.5) 

and (1.2) becomes a pendulum system, for which we know its dynamics. Thus, starting 
with a 2-dimensional subspace II of the codimension 2 manifold X, we are in a position 
to follow the trajectories of (1.2) and study their behavior in the infinite dimensional 
space. 

The phase space of (1.2) can be described in terms of the spectrum of the linear 
operator (for a detailed description see [2], [5]): 

The SG equation arises as the compatibility conditions of the following Lax pair of linear 
operators 

L^(u,C)0 = O, LW(u,C)0 = O (2.7) 

where 

/ is the identity matrix, u = (u(x,t),u t (x,t)) the potential and ( 6 C denotes the 
spectral parameter. 

The spectrum of lS x \ denoted by 

a(L {x) ) := {( e C : L {x) <P = 0, |0| < oo, Va;}, (2.9) 

characterizes the solution of the SG, the potential u also satisfies the SG and is of spatial 
period L. To achieve the spectral analysis of L^ x \ we use Floquet Theory as follows [5]: 
Let us consider the fundamental matrix M(x, x ; u, () of L^> as: 

L^(u,C)M = 0, M(x ,x ;u,C) = / (2.10) 

and the Floquet discriminant 

A(u, C) := traceM(x + L, x ; u, () 
5 



The spectrum of is given by the following condition 

(7(L(*>) := { C G C : A(u,C)eR, |A(u,C)|<2} (2.11) 

where the discriminant A is known to be analytic in both its arguments and is invariant 
along solutions of the SG equation: 

|a(u(0,C) = o 

Thus for each ( A is a sine-Gordon constant of the motion. Moreover, by varying (, we 
generate all of the sine-Gordon constants. Thus, the level sets M u are defined through 
A: 

M u = {vG X : A(u,C) = A(v,C), VCGC} 

We have the following elements of the a(L^ ) which determine the nonlinear mode con- 
tent of solutions of sine-Gordon equations and the dynamical stability of these modes: 

(i) simple periodic (antiperiodic) spectrum (7+ (erf.) 

H A 

<4 = {C 6 C: A(u,C)| c . =c .=±2, -^{,^0} 

(ii) double points of the simple periodic (antiperiodic) spectrum 
^ = {C^ C: A(u,ct = ±2, ^ = 0, ^\^0} 



One studies the stability of a spatially indepedent solution u(x,t) = c(t) of the sine- 
Gordon equation, a solution which is periodic in t. As is well known [7], an analytical 
formula for this solution is 

u(x, t) = c(f) = 2sin~ 1 (msn(t; m)) (2.12) 

where sn(t; m) is the Jacobi elliptic function, and its modulus m measures the amplitude 
of oscillation. Since c(t) is indepedent of x, the linearization can be solved by a Fourier 
transformation from x — ► k = kj = 2nj/L, 

U(x, t) = exp[ik jX } U{kj, t) (2.13) 

This analysis is yields the following result: modes U(kj,t) are exponentially growing in 
t iff fc| e ( 0, m 2 ). That is, the x-indepedent solution c(t) is unstable to long- wavelength 
perturbations. The number K of unstable modes increases with the period L. We list 
the unstable modes 

2tt/L, 2(2tt/L), . . . , K(2n/L) = [m] 
6 



where [m] denotes the largest integer multiple of 2n/L. 

As an example let us compute the spectrum for the spatially and temporally uniform 
solution u = (0,0). In this case, the first system of (2.7) has constant coefficients and 
is readily solved. The Floquet discriminant is given by 



A(u,0 = 2cos [(C + ^t)- 



(2.14) 



Making use of the definition above, the continuous spectrum is given by the entire real 
axis as well as the curve |^| 2 = ^ in the complex £-plane. The simple periodic spectrum 
is given by 



1 (jar 
2V L 



± 



v L 



e z 



(2.15) 



Each of these points is a double point embedded in the continuous spectrum and becomes 
complex if 



0<(^) 2 <1 



(2.16) 



The condition (2.16) is exactly the same as the condition for linearized instability (see 
section 3 below) and the complex double point is given by 



Cd = ^expfi/3] 



(2.17) 



with 



P = tan^ 1 



2jnc 



The Floquet discriminant is an invariant of the SG equation. This means that A(u, C) 
satisfies the following Poisson bracket conditions 

{A(u,C),A(u,C')} = 0, VCC'GC, {A(u,C),^(u)} = 0, VCGC (2.18) 

and the grad of A is given by (see Appendix) 



-(u,C) = ^trace 



M(L)D x [M- l (x)lM(x)] + —M(L)M^(x)£M(x) 

4< 



<5A . . i 
^ (U ' C) = -4^ traCe 



M(L)M- 1 (x)XM(x) 



(2.19) 



1 
-1 



e~" 
s iu 



where M(L) = M(L,x ,u,()- Making use of the fact that (0,0) satisfies the SG equa- 
tion and the spectrum a(L^) is invariant in t, provides a countable number of constants 
of motion. Let u(x,t) denote a solution of the SG equation which is periodic in x with 



period L and U (x, t) denote one of the Backlund transform of u as defined by (A. 2) in 
the Appendix. Then their Floquet discriminant satisfies 

A(C,I/) = A(C,u), VCeC 

that is U and u lie on the same isospectral set in X . When u G X N q (the iV-phase 
quasiperiodic waves), one component of the level set M u is an iV-torus = M u x X NQ . 
Let u G T N have a fully complex double point as an example we return to the case 
u = 0(= (0, 0)) with (d = \e il3 . Let U G M u as defined by 

U (x, t) = 4tan _1 (tan (3 cos (xc cos (3) sech (t sin /?)) (2.20) 

Then [2] 

i. U<£T N 

ii. U is homoclinic to T^, i.e there exists U ± G T N such that || U(x, t) — U ± (x, t) \\ — > 
as t — > ±oo. 

The expression (2.20) is homoclinic to the torus, 

U(x,t) = 0[exp(- v / l-c 2 |t|)], as t — > ±oo (2.21) 

The solution (2.20) is a homoclinic orbit of the integrable SG equation and is related to 
the famous soliton type solution known as the "breather" [2], in the sense that it may 
be called a "spatial breather" , since it is localized in t and periodic in space. 

3 An Invariant Manifold Theorem 

In this section, we state and give the proof of an invariant manifold theorem for the 
perturbed sine-Gordon equation. There exists a two dimensional subspace X 

X = [ueM 1 : ^ = 0} 

such that the set II = X fl H 1 is invariant under the flow of the perturbed system. The 
manifold II is real symplectic with a nondegenerate 2-form flu which arises from the 
restriction of the symplectic structure to II. For e — 0, the system (1.2) restricted to 
II becomes a 1 degree-of-freedom completely integrable Hamiltonian system, pendulum 
system, for which we know its dynamics. Furthermore, the uniform solution O = (0, 0) 
is a saddle type for the reduced pendulum system which lies on the torus S. On the 
whole space this uniform solution becomes a singular point of type saddle-focus. The 
breather solution U (x, t) of the full system approaches this set when t — > ±00. Starting 
with the space II, we are in a position to follow the trajectories of (1.1a) and study the 
behavior of its solution in the infinite dimensional space. 
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3.1 Linearized Analysis 



The linearized system of the perturbed sine-Gordon in the neighborhood of the solution 
O — (0, 0) e S takes the following form for the Case I; 

r t = Ar + Kj(r) (3.1) 

where r = (f , f t ) 

A = ( c 2 d 2 + 1 ) ' K/(r) = ( -eaf t - G T (r; e) ) 
and the function Gi is defined as follows: 

G/(f) = egi (O + f) + sinf - f = 0(r 3 ) (3.2) 

for the Case II, 

r t = Ar + K n (r) (3.3) 

where 

Kn ^ ( -eaf t - e) ) 

and the function G/j is defined as follows: 

G n = eg 2 {t) + sinf -r = 0(f 3 ) (3.4) 

We can also combine the two equations (3.1) in a single equation, (respectively for the 
Case II) 

ftt ~ c 2 f xx - f + eaf t = Gi(r) (3.5) 

Let e — 0, (3.5) becomes 

f t t - c 2 ^ - f = 0, c 2 e(i,l) (3.6) 

Substituting f = f (t)e lfcjX with fcj = and j an arbitrary integer we obtain from (3.6) 

f" - (1 - c 2 k])f = 

This shows that the j-th mode grows exponentially if c 2 {2nj / L) 2 < 1, j — 0, 1. Using 
now the even spatial periodic boundary condition we see that our problem has two 
unstable directions, corresponding to the unstable directions of the saddle point (0,0): 
The zeroth mode in the plane n a — 1, which is the most unstable in the sense that it 
grows at the fastest rate and j — 1, «i = — {2nc/L) 2 off the subspace II. Hence, the 
saddle point O = (0,0) has a 2-dimensional unstable manifold. If c 2 (2irj/L) 2 > 1, we 
find infinite number imaginary eigenvalues Xj = iictj = ±i^/ c 2 (2nj/L) 2 — 1 , j > 2. 



In order to study the local behavior of solutions of the nonlinear system (3.1) we 
need to analyze the spectrum of the operator 



Ce ~ ( c 2 dl + 1 -ea ) 

We consider the eigenvalue problem C e e = Ae for the eigen-pairs {e(x),X}. Using 
Fourier expansions of e(x) in exp(ikjx) leads to a quadratic expression for the eigenvalues 
A: 

A 2 + \ea - (1 - c 2 kj) = 0, kj = 2vrj/L, j = 0, 1, 2, . . . 
For j = 0, 1, we have: 

e s ,u = (l^TCTj^coskjX 



while for j > 2, the eigenvalues come in complex conjugate pairs. 
3.2 Invariant Manifolds-Case I 

From the above analysis, we conclude that the function u may be written as a linear 
combination 

u = r u e u + f s e s + f c 
and in terms of these variables the equation f t = C e r splits into 

f Cjt = Af c (3.7) 

where the operator A corresponds to the infinite number of center directions. In a 
neighborhood of the fixed point (0, 0) which lies on a torus S, the nonlinear equation 
(3.1) can be viewed as a perturbation of the linear equation (3.7). Under the flow of 
this linear equation and for e = 0, S has a 2-dimensional stable and unstable man- 
ifold together with a codimension 4 center manifold. We focus our attention on the 
center manifold E C (S) together with center-stable E CS (S) and center-unstable E cn (S) 
manifolds, on which the solutions have growth rates bounded by 

(i) exp[(jjt/n], for t > 0, and exp[—<jjt/n], for t < 0, j = 0, 1 



(ii) exp[<7j|i|/n], for Vt, j > 2, aj = iy (c/cj) 2 — 1 
respectively, while for the operator A we have 

||exp[^]|| <nCexp[aj\t\/n] 
10 



for some integer n and positive constant C. 

Let us now concentrate on the existence of locally invariant manifolds of the per- 
turbed SG equation. To this end, we introduce a localization function, [16], such that for 
all Ci, c 2 with < C! < c 2 and integer p > there is a function ip§ G C p , (ps(y) = ^{y/S): 



if : R > R, <p(y) = 





for 


\y\ 


< 


Cl 


{! 


for 


\y\ 


> 


c 2 



with supJ</?W(y)| < c for all 2 < p and suitable c. Hence, equation (1.1a) is written as 
a perturbation of (3.7) with u = (f s ,f u ,f c ) 

r u ,t = X £ u f u + P u (f,5;e) 
r s ,t = -Xlf s + P s {r,5;e) 

f c ,t = Ah + P c (f,S;e) (3.8) 

where P u , P s , P c are nonlinear functions. This localisation has the effect of keeping the 
flow unchanged in the neighborhood of torus S. We are now in position to prove the 
existence of an infinite dimensional, locally invariant manifold for equation (1.1a) in the 
vicinity of torus S, local invariance means that solutions can only leave the manifold 
through its boundary. The manifold we identify for equation (1.1a) is a codimension 4 
center manifold for the torus S for e — 0, and continuous to exist for nonzero values, even 
through S is usually destroyed by the perturbations. We also show that the manifold 
possesses locally invariant stable and unstable manifolds which codimension 2. 

The following theorem concerns the existence of the locally invariant manifolds de- 
scribed above. 

Theorem 3.1 There exists a 5 neighborhood of S C H 1 and eo(S) > such that 
\/e G [0, £o) equation (3.8) has locally invariant center-stable and unstable manifolds of 
codimension 2 

W? oc {M e ) = {rem 1 : r u = (f s ,f c ;e)} 

W? oc {M e ) = {fGtf:r; = h s (f u ,f c ;e)} (3.9) 

with f u = (r U) o, r U; i), f s = (f S) o,f Bj i) and a codimension 4 center manifold 

M £ = {reM 1 :f u = h cu (r c ;e),r s = h cs (r c ;e)} (3.10) 

where f c is defined all ofR and satisfies sup teR ||r c || H i < while 5' belongs to a neigh- 
borhood of 5. 

Proof. We write (3.8) in the form of integral equations 

f u (t) = exp[\ £ u (t-t u )]r u (t u )+ f exp[A £ u (t-0] J P u (r(0,5,e)de 

Jt u 
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f s (t) = exp[-Xt(t-t s )]r s (t s )+ [* eM-K(t-Z)]Ps(m,6,e)dt 

f c (t) = exp[^]f c (0) + [* exv[A(t-Z)]P c (r(Z),6,e)dt (3.11) 

Jo 

From the gap in the growth rates and the definition of the flow F*, we characterize the 
invariant manifolds Wf oc (M £ ), W^ C (M £ ) by 

W* oc {M £ ) = {reM 1 : sup t > {exp[- C r 1 t/ K ]||F*(f;£)|| Hl }<oo} 
W? oc (M £ ) = jf G H 1 : sup t < { exp[o"it//t ] ||F*(f; e) || H i} < ooj 

Focusing our attention upon W^ C (J\4 £ ), for f in a sphere B of arbitrary radius g, we 
introduce the norm 

\\r\\n : = sup{exp[-<7it///]||f|| Hl }, g < \\f\\ ml 
t>o 

For f G Wf* c (M £ ), we have exp[— A^ u ]|f u (i u )| — > as t u — > oo. Thus, on the center- 
stable manifold the integral equations (3.11) can be written as: 

fu(t) = f exp[A£(f-0]Pu(f(0,*,e)d£ 

f B (*) = exp[-A^]f s (g+ /' exp[-A|(*-0]P B (r(0^,e)de 

•/ o 

f c (t) = exp[^]f c + /%xp[^(*-0]P c (r(0,^)d£ (3.12) 
To prove the existence of Wf oc {M. £ ), we use Newton's iterations. Let f° = and 
C +1 (t) = /" exp[A^-0]P u (r m (0,Md£ 

J +oo 

C +1 (t) = exp[ - Kt\f s + /* exp[ - K{t - 0] P s (r m (0, 5; 5) d£ 

f c m+1 W = exp[^]f c + /* exp[A(t-0]Pc(f m (0^;e)d^ (3.13) 

Jo 

We note that P(f, 5; e) is a smooth function whose terms are linear with coefficient 
e or nonlinear and localised in a (^-neighborhood of 5. If we let P' be the derivative of 
P, we obtain the following estimate 

\\P(r,5;e)\\ m <\\P'\\ \\f\\ m +e 

where \\P'\\ is the supremum of the magnitude of P' which is equal to C(S + e). 
For if ||r m || Ko < C, we obtain from (3.13) 

l|f m+1 (OII]Hi < «oC , exp[<Ti*/K ](||f s || H i + ||r c || H i + e) 

12 



+ f exp[ y {t " M5 + £) Ir(0 IIhi de 

+ /* k C(5 + e)exp[ - HIHOIU ^ 

By using the bound of f m we obtain 

||r m+1 (t)|| el < C[||f s || Hl + ||f c || Hl +£^(5 + £)||p 1 || Ko ]exp[a 1 t// t0 ] 

where the constant C is independent of K ,e and <5. Fix 5 — where a = C/4, we 
obtain for all £ < C/4kq: 



r m+1 (*)IL <C(p) + -||f 



m II 
2 ii- ll*eo 



Thus, we have proved that the sequence r m is well defined and ||r m || Ko < 2C(p). 
Since, the nonlinear term is smooth, we have a similar estimate for the difference 
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|~m+l_~m|| < -||f m — f m_1 || 



(3.14) 



which implies that r m — > r and 

||r|| K0 < 2C(£ + ||f c || Hl + ||f s || Hl ) 

We note that all terms in (3.12) are smooth, which entails that {f m } is differentiable. 
We observe that the derivative Df m satisfies: 



exp[y(t-0 
+ C f exp[^(t -0}\\P'Dn\m^^ 

JO l ^Ko J 

By using the bounds on P', we obtain 



\P'Dr T 



Ih 1 



\Dr 



m+l| 



< C + -||£>f m || 

ko — 2 



which implies ||7Jr m || < 2C. 

r II II /t — 

Now, we are in position to estimate the difference between two terms in the sequence 
{f m }, 5f m = f m — f m ~ l . We note by the mean value theorem 



p\n - p'(f-y\\ ml < cw 



Ih 1 ll r He 1 



We have: 



\\D5f 



E 1 



< C f 'exp[^(t-C)}{\\P'(D5n 



Ib 1 
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+ {\\i^\\^ + \\i^- 1 \\ m )\\s^\\ m }^ 

+ cj\ w [^-(t-0]{\\P'(D5n\^ 



(3.15) 



and 



Then, we obtain 



\\6n\mA\Dn 



E 1 



< exp 



2ai 



l^mn II £^11 

I 1 1 K() 1 1 II KQ 



£>f m+1 - Df E 



Iko/2 



< C||(Jf 



+ - \\Df m — Dr m ~ l \ 
k 2 11 1 



and sequence {r m } converges in C 1 using the ||-|| Ko / 2 -norm. We repeat this procedure 
to find bounds on {D j f m } in the ||-|| K ^.-norm with n/j > 2. 



We obtain r G C r for r < 



— 1. Therefore, we define 



f u = h»(f s , f c ; e) = f exp[ -A u £ ]^u(r(0, 5; e) d£ 

which is a C r functional with || -Df u || < 1/2 and 

WL(M £ ) = { f G H 1 : f u = K 1 ( r~ s , f c ; e) } 

is a C r manifold of codimension 2, with f u = (f Uj o, ^ u ,i)- Exactly analogous results are 
obtained for the center-unstable manifold, W^ c (Ai £ ). 

The center manifold M. £ is described by W^.Me) H H^" C (A^ £ ), and equivalent we 
need to solve the following system 



r u = h u (r s ,r c ;e), r s = h s (r u ,r c ; e) 



(3.16) 



we note that ||D/i u ' s ||<l/2 and (3.14). From the implicit function theorem, we obtain 
that the system (3.16) has a unique C T solution given by 

f u = h™(f c ;e), f s = hf(f c ;e) 

□ 



3.3 Invariant Manifolds-Case II 

The time-periodic perturbed sine-Gordon equation (1.1b) can be rewritten as: 

U( = Au + B(u) + e(Cu + g 2 (t)) 



(3.17) 



14 



with 

A:=f ° M, B(„):=( ° V C : = f ° ° 
V c z oz \ smu / \ —a 



X 



We assume that g 2 : H 1 x S 1 -> H 1 is C 00 T-periodic in time where S 1 = R/T. The 
associated extended autonomous system in a cylinder in the space H 1 x S 1 

ii = Au + B(u) + e{Cu + g 2 (6)) 

= 1 (3.18) 

has a smooth flow F*. The geometry of the unperturbed system (3.17) is the same with 
the case I. 

Suppose that linearized system of (3.17) with respect to the singular point p = O 

r t = C £ v + eg 2 {t) (3.19) 

with 

C - ( ° 1 ) 

^ - { c 2 d 2 x + 1 -ea J 

has a T-periodic function f(t,e), such that f(t,e) = 0(e). 

For e > 0, the operator e TCe has a spectrum consisting in two pair of real eigenvalues 
for j — 0,1 and infinite number of complex eigenvalues for j > 2, spec(e T££ ).,>2, in such 
a way: 

Hie < dist^spec(e T£e )j> 2 , \z\ = ij < [i 2 e (3.20) 

with fj,i, /i 2 > 0. 

We denote by P e : H 1 — > H 1 the Poincare map for the flow F* such that P°(po) = Po 
and p £ the fixed points of P e correspond to periodic orbits of F*, near to p = O : p £ = 
Po + O(e)- 

Indeed, by (3.19) we have: 

r(t, e) = e tCs r(0, e) + f ee^ )Ce g 2 (0 d£ (3.21) 

Jo 

and r(T,e) = f(0, e), \\r(t, e) ||<//3£. 

We seek a solution u(t,e) of the system (1.1b) such that: 

u(t,e) = e^u(0,e)+ f ee^^g^) d£ + /%e(*-^(B(«(£,e)) - u(£,e))d£ (3.22) 

JO JO 

From (3.21) and (3.22) we obtain: 

u(t, e) - f(t, e) = e^ (u(0, e) - f (0, e)) + f ee^^(B(u^, e)) - u(£, e)) d£ (3.23) 

JO 

Thus, 

\\u(t, e) - f(t, e) || Hl <^ 4 £ + / IKf , e) IU d £ 

«/ 
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using \\f(t, s)\\ m i<(J>3£- By Gronwall's estimate we obtain \\u(t, e)\\ ml < /iqE. 
We consider now the set B £ as 

B e = {f{t,e):\\f(t,e)-f{0,e)\\ Kl <e} 

and the map 

P £ : B £ -> H 1 
u(0,e) -> u(T,e) 

We seek a fixed point of P £ . From (3.23) u(0, e) is a fixed point of P e if and only if it is 
a fixed point of the following map 

G e : B £ -> H 1 

with 

G e («(0, e)) = f(0, e) + L7 1 jT e^ (b(u(£ e) - u(£, e))) d£ (3.24) 

where L £ = (Id — e T,C£ ) is invertible operator and the condition 

ULrini <ine~ l (3.25) 

is inferred from (3.20). The conditions (3.24) and (3.25) entail 

||G e («(0,e))-f(0,e)|| Hl <e 

which implies that for e <^ 1 G £ maps B £ to itself. Also, G e is a contraction since, for 
e small enough 

||L>G £ K0, £ ))|| Hl = || L; 1 [ T e( T -^D(B(u(Z,e)-u(Z,e)))oDFl(u(0,e))dZ\\ 

JO s H 1 

< H%E < 1 

Thus, G e has a unique fixed point on B £ . There is an analogue result for the Poinacre 
map Pf with the section H 1 x {t } in H 1 x S 1 . 

In the following Lemma, we describe the geometry of the invariant manifolds for 
the non-autonomous perturbed sine-Gordon equation (1.1b). Its proof based on the 
invariant manifolds theory [12] and the smoothness of the flow operator F*. 

Lemma 3.1 Let p £ (t ) denote the unique fixed point of Pf Q . Corresponding to the 
spectrum of e tCs , there are unique invariant manifolds W cs (p e (to)) (center-stable) and 
W n (p £ (to)) (unstable manifold) for the map Pf such that they are invariant under the 
Pj o? tangent to the eigenspaces of e tCe respectively at p £ (t ) and C r -close as e — > 
to the homoclinic orbit U(x,t) of the sine-Gordon equation (1.2) for t s ,t u such that 
t s < t < 00, —00 < t < t u , respectively. We denote by q e (t) = (p £ (t),t),t G (0,T) the 
periodic orbit of the extended system (3.18) with q e (0) = (p £ , 0). The invariant manifolds 
for the periodic orbit q £ (t) are denoted W cs (q £ (t)), W u (q £ (t)) and 

W j (p 6 (t )) = W j (q £ (t)) n (H 1 x {t }), j = cs,u 
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4 Mel'nikov Analysis 



In this section, we shall combine the invariant manifold theory of the previous section 
with explicit global information from the unperturbed sine-Gordon equation to establish 
the existence of persistent homoclinic orbit to the singular point £ . 

4.1 Geometry of Mel'nikov Method-Case I 

We recall that the unstable manifold W u (0 £ ) C W n (A4 £ ) is two-dimensional with one 
direction in the plane of spatial-independent solutions II and the other direction off the 
plane in the function space H 1 . We consider an orbit in W n (0 £ ) that leaves £ and 
takes off away from the plane, then defines an orbit which lies on the codimension 2 
center-unstable manifold. We will construct a distance function T> between the stable 
and unstable manifolds, whose zeros correspond to orbits that do not lie on the plane 
II and are asymptotic to the saddle-focus point £ in backward and forward time. 
Let us rewrite equation (1.1a) in the vector form 

u t = JV u #(u) + eGj(u,fi) (4.1) 

where G/( u, /x ) = (0, bgi(u) — au t )) ± and /i = (a, b, c) 

For e < eo, (3.8) is related to the original system (1.1a) within some fixed open set 

Y = { (r s ,f u ,f c ) G H 1 : |f s | < C s , |f u | < C u , ||f c || Hl < C c ) 

where C s ,C n , C c are fixed positive constants. We are interested in the dynamics of 
the solution u £ {t) = (r s ,f u ,f c ) of (3.8) in forward and backward time. We define a 
neighborhood U of M £ 

U = { (f s ,f u ,f c ) G Y : |f s | < 5, |f u ] < 5, \\f c \\ m < 5 < C fc } 

and its boundary dU consists of two parts dU s stable and dU u unstable as follows 

dU s = {(r s ,f u ,f c ) eU : |f s | = 5} 

dU u = { (f s , f u , f c ) eU : |f u | = 8 } 

where the solution curves enter through dU s and they exit through dU u . As was said 
before, the breather solution U (x, t) is also a solution of the SG and is homoclinic to the 
point O = (0, 0) G II which has two unstable directions, (cf. section 3). The unperturbed 
system possesses a homoclinic orbit u which evolves along the unstable direction from 
the plane II, intersects the unstable boundary dU u at the point Uq and defines a unique 
orbit (u n (t),t) which lies on the codimension two center-unstable manifold. After a 
finite time, this homoclinic orbit u intersects the stable boundary dU s at the point 
«q and comes back to the plane II C S and is asymptotic to the singular point O, 
following a unique curve (u s (t),t) which lies on the center-stable manifold. Thus, for 



17 



e 7^ 0, there exists a perturbed orbit u £ which evolves from the perturbed saddle point 
£ = O + 0(e) E II and intersects the unstable boundary at the point u^. 

Now, choose a point Uh on the homoclinic orbit for the unperturbed system. Let 
E be a codimension 2 hyperplane at Uh which is transversal to the homoclinic orbit u 
and which contains the vector V = (^, JA). W s (M.q) fl E has a codimension 1 in E 
and z7 is transversal to W S (M ) fl E. Let w s be the intersection of the line through 
u u e W u (0 £ ) C W^ u (7W e ) and the manifolds W S (M £ ) along to direction V. We may 
thus define the inner product 

V= (WA,u u -u s ), (4.2) 

as a measure of the distance between u u and u s . 

To calculate D, we parametrise the orbits as follows, 

for t<0, u = u h (t + Q, u n = u*(t + Q 

for t>0, u = u h (t + Q, u s = u s £ (t + t*) 

where u s ,u u be the solutions of the localized equations (3.8). Since Uh remains in a 
^-neighborhood of S for t > t* u is a solution of the localized system (3.8), we recall 
that the invariant manifold theorem holds in a ^-neighborhood of S C H 1 . For t< — t* 
both orbits u, u u remain in a 5-neighborhood of S. We apply the Gronwall's estimate 
to system (3.8) and from ||w u — Uh\\ Ml < C\E, for t < we have, 



' u h \\ m < eC 2 e- at (4.3) 



Similar, for t > 



u 



u h \\ ml < eC 2 e° l (4.4) 



The above analysis allows to decompose (4.2) as follows; 

V = (VA(u h ),u u - u h ) - (VA(u h ),u s - u h ). 

Let 

V s = (X7A,ul(t) -u(t)), t>0 
V u = (VA,<(t) -«(*)), t<0 

Then 

D = V u (0) -V s (0) 
Proposition 4.1 The distance function (4-2) is given by 

V = V U (0) -V s (0) =eM I + 0(e 2 ) (4.5) 

where 

M I= / (VA(«(t)),Gj(u(*)))dt, (4.6) 
denote the Mel'nikov integral with (a,b) = Jq a(x)b(x) dx. 
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Proof. We start the proof by asserting from (4.1) that u = u(t), u £ = u £ (t), u s £ = u s £ (t) 
solve the following system of equations 

d t u = JVH(u) 

d t u u £ = JVif(<)+eGj«,^) 

d t u\ = JV^(<)+eG/«,Ai) (4.7) 
Differentiating V u (t), we have 

Z>«(t) = e (VA(u), Gj(u)) + (V 2 A(u) JVif («),<- u) 

+ ( VA(«), V(J Vif («)) « - u) ) + (VA(tt), 5«, u) ) 

where 

J Vif «) = J VH(u) + V( J Vff(u)) « - u) + S(u u £ , u) 

we leave out the argument t of the functions for the sake of a less cumbersome notation. 
From (2.17) we have (JVH, VA) = 0, differentiating with respect to u gives 

V(JV#) X VA + V(VA) JVH = (4.8) 

where - 1 denotes the matrix transpose. Taking the inner product of (4.8) with — u 
gives 

(V 2 A J VH, — u) + ( VA, V(JVif)(<-u)) =0 

and thus, 

X> u (t) = e{ VA(m), Gj(«) ) + ( VA(u), 5«, «) ) (4.9) 
There exists > for t G (— oo,i*) such that 

||VA(u(f))||Hi<Ce rt , a:=y/l- (4.10) 

since the VA is evaluated on the homoclinic solution which satisfies the condition (2.21) 
and (4.3) entails 

||«S|| H1 < IK - «||hi < Cie- 2 V (4.11) 

with 25 < a. Now, we want to find out the limit \im t ^_ OQ V n (t) in order to represent 
V(0)hy 

X? u (-oo) + [° V u (t)dt (4.12) 



We calculate lim t ^_ 00 P u (t) and show that lim^-oo T> u (t) = 0. We have 

lim V u (t) = lim (VA(u(t)), <(*) -«(*)) 

t^ — OQ t— » — OO 

There are C 2 , C 3 > such that 

K(T-%)-KM\\m < c 2 ^ r 
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\\u(r-e )-Fl(u)\\ ml < C 3 e A - (4.13) 

for all r G (-00, 0], where Ai > 0, u^(-t^) = v%, u(-t%) = u u and 

F T £ (u) ^0 £ as r -> -00 (4.14) 

FJ(«) = O G S* (4.15) 
By these relations (4.13-4.15), we obtain that as £ — > —00 

K(t)-ii(t)llm~0(f) (4.16) 
From relations (4.10) and (4.16) we see that 

lim £> u (t) = 

All this finally yields upon integration of (4.9) 

V u (0) = e f (VA(u(*)), Gj(«(*))) dt + 0(e 2 ) (4.17) 

Similarly, we have 

POO 

V s (0) = -e / (VA(u(t)), Gj(u(t))> dt + 0(^ 2 ) (4.18) 
Jo 

whence the distance function (4.5) is given by 

D = D u (0)-D s (0) = eM/ + O(£ 2 ) (4.19) 

with ,00 

M I =\ (VA(u(t)),Gj(u(t)))dt, f = 1,2 (4.20) 

j —00 

where w(t) = U(x,t) the homoclinic solution of the unperturbed system. 

The Mel'nikov integral is convergent. Because, the perturbation term Gj = bgi(u) — 
au t is exponential decreasing expression, since we have assumed that g\{u) satisfies the 
conditions <7i(0) = 0, <9 u gi(0) = and in the neighborhood of zero 

lkiHII H i</i9||«|| 2 </i 9 e- 2<J * (4.21) 

evaluated on the homoclinic orbit, with fj,g > 0. Also, hold 

\\u t \\ m i<l^ioe~ at 

and (4.10). Thus, obtaining convergent integral M. □ 

Using the integral of motion, we may derive for the case I, explicit formula for the 
Mel'nikov function 

/oo 
(VA(U(x,t)),Gj(U(x,t))}dt (4.22) 
-00 
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evaluated on the homoclinic orbit U(x,t). Thus, 

roc rL §/± 



roc rL A/\ rOO rL fi^ 

M I (p,n)=b / —(U)gAU)dxdt-a / — (U)U t dxdt (4.23) 

J-oo JO OUt J-oo JO OUf 



with (5 defined in (2.17) and 

5A 



i ( 1 



5u t Ac\2$m2f3 



(4.24) 



where E, H are smooth functions of the general solutions ipi,ip2 of the Lax pair (2.7). 
Setting Mj( /?,//) = 0, we obtain an algebraic equation for parameters 



b — an = 



(4.25) 



where 



roo rL / roo rL A /± 

We denote the surface defined by (4.25) by £p 

£p : P = B(b ,a ;c) 
There exists a region TZ of the surface £p such that 



-L 5 A 



d(3 



(Po,bo,a ;c) ^ 



(4.26) 



and 



df3 



< I, for I > 0. By the implicit function theorem there is a neighborhood TZ 



of (bo, a ; c) and a unique C n 2 function 

Po = B(b , a ; e, c) 

defined in TZ such that 

B(b Q , a ; 0, c) = p 

and V(B(bo, ao; £, c), &o, «o; £, c) = 0. Since "D (the distance of stable and unstable man- 
ifolds) is (7 n ~ 2 smooth function by (4.26), we obtain 

d 



dP 



■V(B(bo, a ; e, c), b , a ; e, c) ^ 



and 



df3o 



< m, m > 0, for (& , ao; c) G TZ. Then, W s (yVf £ ) and W U (.M £ ) intersect 
transversely in the neighborhood of TZ. 

Proposition 4.2 There exist a neighborhood TZ for the parameter b, a, c and a positive 
constant eo > such that fprany \e\ < eo and (b,a,c) G TZ the W s (J\4 e ) and W n (M. £ ) 
intersect transversely. 
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4.2 Case II: Time-Periodic Perturbation 



The system (1.1b) can be rewritten in the more general form as follows: 

x = f(x)+eGjj(x,0), 

9 = 1 (4.27) 

We assume that G n G (^(H 1 x S^H 1 ), where S 1 = R/T. This implies that the 
associated suspended autonomous system on H 1 x S 1 has a smooth local flow F* : 
I 1 x S 1 -> H 1 x S 1 , this means that F* is a smooth map defined for small \t\ which is 
jointly continuous in all variables e, t G R, x G H 1 , 9 G S 1 . 

We recall that for e > small, there is a unique fixed point p £ of P £ near p = O = 
(0, 0), moreover p £ — po = 0(e), i.e there is a constant K such that ||p e || < Ke, Ve > 0. 

We consider a neighborhood of the codimension 2 center manifold Ai £ , in which we 
attempt to establish the transversal intersection of W s (p £ ) C W S (M £ ) and W n (p £ ) C 
W U (A4 £ ). Based on the assumptions and preliminary results of Lemma 3.1, we pro- 
ceed to calculate the distance T> between the perturbed manifolds W s (p £ ),W n (p £ ), by 
calculating the 0(e) components of perturbed solution curves of equation (4.27) from 
the first variation equation. We expand solution curves in W s (q £ ) and W n (q £ ), as we 
described in Lemma 3.1, points in W s ' u (p £ ) are obtained by intersecting W s ' u (q £ ) with 
the section H 1 x {t }. 

Therefore, (4.27) possesses the perturbed orbit 

x^, t ) = ^(t - t ) + exy(t, t ) + 0(e 2 ) 

on the perturbed unstable manifold, where x"(t, t ) is the solution of the first variation 
equation 

x?(t, t ) = Di(x (t - * )K(*, to) + G//(x (t - t ), t) 

and a similar one lying on the perturbed stable manifold. 
As in the case I, the distance function is given by: 

V(t ) = eM n (t ) + 0(e 2 ) (4.28) 

with 

M n (t ) = r [ L ^(U) (g 2 (t - t ) - aU t ) dxdt (4.29) 

J~ooJ0 OUt 

evaluated on the homoclinic orbit U. 

For g 2 (t) = Tcosut, the equation (4.29) takes the form: 

roo r-L A A 

M n (t )= / — (U) (Fco8u(t - 1 ) -aUtjdxdt (4.30) 

J-ooJO 0U t 

The Mel'nikov integral in equation (4.30) is convergent. We therefore consider the 
improper integral as the following limit of a sequence in time {r^' u } = {±2nn/u}, 
n = l,2,... 

,,2™ L ^ 

M n (t , T, a, 0, c) = Jim J ^ ^ J^( U "> ( ^^(t - *o) ~ aU t ) dx dt (4.31) 
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Since jj^(U) and U t are rapidly decreasing expressions, we are allowed to extend the 
integration limits in (4.31) to infinity and thus obtaining convergent integral. 
We find that the Mel'nikov function is of the form 

M n (t , T, a, (3, c) = ^ 2 rsinwt - ^2 (4.32) 

where 

r°° r L 5 A 

A2sincuto = / - — (U(x, t))cosu(t — to) dx dt 

J -oo Jo ou t 

1 2 = — (U(x,t))U t (x,t)dxdt 

J -oo Jo ou t 

Setting Ma (to) = 0, we obtain 

t^±-sm-\^ X 2) X2 = ^- (4.33) 

LO 1 J\2 

Lemma 4.1 For | < c 2 < 1, there exists an open interval I C K, such that for any 
T,a,(3 G /, there are two values of t : to = ^o( c ? u i T, a, (3, X2) (cf (4-33)) at which 
Mn(to,n) — 0, /I— (c, uj, T, a, (3). Moreover, when M n (t ) is viewed as function of t^ , 
these zeros are simple. 

By this lemma, eq. (4.28) and the implicit function theorem, we have the following 
proposition: 

Proposition 4.3 There exists a subregion TZi of external parameter space 

n= {(c,u,T,a,(3) : - < c 2 < 1, uj G (0,0.87), T,a G (0.1, 1), (3 G ( -7r/2,7r/2)} 

and a positive number So, such that for any fixed parameters {c, uj, T, a, (3} G 7Z± x 
(0,£q); there are t^ (cf. (4-33)) at which W s (p £ ), W u (p £ ) intersect transversely. 



5 Conclusions 

In this paper, we have used a Mel'nikov type analysis to detect the splitting of homoclinic 
manifolds in a perturbed sine-Gordon equation with even spatial periodic boundary 
conditions. We extended the two-dimensional phase portrait for the pendulum into the 
infinite dimensional phase space H 1 of the integrable periodic sine-Gordon equation. 

This analysis consists first of a study of homoclinic and chaotic dynamics on per- 
turbed SG systems. The prerequisite for such study is that the unperturbed soliton sys- 
tem, viewed as a Hamiltonian system should have a homoclinic structure. In the process 
of proving the existence of invariant manifolds for such a system several mathematical 
tools from dynamical systems theory are utilised. We established, via Mel'nikov theory 
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the existence of transversal intersection of a codimension 2 center-stable and center- 
unstable manifold, of the perturbed singular point O e , which implies the occurrence of 
chaotic dynamics. In particular geometric settings, a simple zero of the Mel'nikov func- 
tion, with respect to one of its parameters insure the intersection of stable and unstable 
manifolds of a critical torus and the persistence which follows as a consequence of this 
intersection. 

Acknowledgement 

The author had stimulating discussions with T.Bountis and G.Haller. He would also like 
to thank the anonymous referee for helpful comments and remarks. This research was 
supported by the postdoctoral fellowship No. ERBFMBICT983236 of the Commission 
of the European Communities. 

References 

[1] Bishop, A.R., Forest, M.G., McLaughlin, D.W. & Overman E.A.II, Physica D 23, 
pp. 293-328, 1986. 

[2] Ercolani, N.M., Forest, M.G. & McLaughlin, D.W., Physica D 43, pp. 349, 1990. 

[3] Ercolani, N.M. & McLaughlin, D.W., Towards a Topological Classification of In- 
tegrate PDE's in The Geometry of Hamiltonian Systems, Springer- Verlag, New 
York, 1991. 

[4] Faddeev, L.D. & Takhtajan, L.A., Hamiltonian Methods in the Theory of Solitons, 
Springer- Verlag, Berlin, 1987. 

[5] Forest, M.G. & McLaughlin, D.W., J. Math. Phys., 23, pp. 1248-1277, 1982. 

[6] Forest, M.G. & McLaughlin, D.W., Stud. Appl. Math., 68, pp. 11-59, 1983. 

[7] Guckenheimer, J. & Holmes, P.J., Nonlinear Oscillations, Dynamical Systems and 
Bifurcations of Vector Fields, Springer- Verlag, Berlin, (2nd ed.) 1986. 

[8] Haller, G. & Wiggins, St. Multi-pulse jumping orbits and homoclinic trees in a 
modal truncation of the damped-forced nonlinear Schrddinger equation Physica D 
85, pp. 311-347, 1995. 

[9] Haller, G. Multi- dimensional Homoclinic Jumping and the Discretized NLS Equa- 
tion, Comm. Math. Phys. 193 (1), pp. 1-46, 1998. 

[10] Haller, G. Homoclinic Jumping in the Perturbed Nonlinear Schrddinger Equation, 
Comm. Pure Appl. Math. 52 (1), pp. 1-47, 1999. 

24 



Haller, G. Chaos Near Resonances, to appear in Springer Verlag, 1999. 

Hirsch, M., Pugh, C. & Shub M., Invariant Manifolds, Springer Lectures Notes in 
Mathematics, 583, 1977. 

Holmes, P.J. & Marsden, J. A Partial Dif ferential Equation with Infinitely Many 
Periodic Orbits and Chaotic Oscillators of a Forced Beam, Archive for Rational 
Mechanics and Analysis, 76, pp. 135-165, 1981. 

Kovacic, G. & Wiggins, St., Orbits Homoclinic to Resonances, with an Application 
to Chaos in Damped and Forced sine-Gordon Equation Physica D 57, pp. 185-225, 
1992. 

Li, Y, McLaughlin, D.W., Shatah, J., & Wiggins, S., Persistence of Homoclinic 
Orbits in Perturtbed NLS Pde's, Comm. Pure and Appl. Math., Vol. XLIX, pp. 
1175-1255, 1996. 

Marsden, J. & Cracken, M., The H op f Bifurcation and its Applications , Appl. Math. 
Sci., 19, Springer, Berlin, 1976. 

McLaughlin, D.W. & Shatah, J., Homoclinic orbits for PDE's in Recent advances 
in partial differential equations, Venice 1996, pp. 281-299, Proc. Sympos. Appl. 
Math., 54, Amer. Math. Soc, Providence, RI, 1998. 

McLaughlin, D.W., Whiskered Tori and Chaotic Behavior in Nonlinear Waves, in 
Proceedings of the International Congress of Mathematicians, Vol. 1, 2 (Zurich, 
1994), PP- 1484-1493, Birkhauser, Basel, Basel, 1995. 

McLaughlin, D.W. & Overman E.A.II, Surveys in Applied Mathematics, 1, 1995. 

Mel'nikov, V.K., On the Stability of the Center for Time Periodic Oscillations, 
Trans. Moscow Math., 12, pp. 1-57, 1963. 

Pedersen, N.F., Advances in superconductivity, Deaver and Ruvolds, eds. NY 
Plenum, 1982. 

Rothos, V.M. & Bountis, T.C., The Dynamics of Coupled Perturbed Discretized 
NLS Equations, Physica D113, pp. 326-331, 1998. 

Rothos, V.M., MeVnikov Theory of Coupled Perturbed Discretized NLS Equations, 
to appear in Chaos, Solitons and Fractals, 1999. 

Temam, R., Infinite Dimensional Dynamical Systems in Mechanics and Physics 
(2nd ed.), Springer-Verlag, NY, 1996. 

Wiggins, St., Global Bifurcations and Chaos: Analytical Methods, Springer-Verlag, 
1988. 



25 



